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1 Introduction 

It is well known that harmonic functions (i.e., solutions of the Laplace equation Au = 0, 
where A = Yli=i Jr) P ossess the mean value property: Namely, if u is harmonic on a 

i 

domain Q C M. d , then for every closed ball B r (x) C Q of a center x G f2 and radius r > 
the average of u over B r (x) equals to the value of x, i.e., 

u(y)dy = u(x), 

B r (x) 



B r X) 



where |5 r (rr)| is the volume of the ball B r (x). There is a similar statement for mean 
values on spheres. More generally, if u is not assumed harmonic but Au is a continuous 



function, then 



\ I / u(y)aly - u(x)) = c n Au(x) (1.1) 
r^or 2 \\B r (x)\ J Br[x) J 

for the appropriate dimensional constant c n . 
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What are the fractal analogs of these results? The analytic theory on p.c.f. fractals 
was developed by Kigami [3j HJ [5] following the work of several probabilists who con- 
structed stochastic processes analogous to Brownian motion, thus obtaining a Laplacian 
indirectly as the generator of the process. See the book of Barlow [I] for an account 
of this development. Since analysis on fractals has been made possible by the analytic 
definition of Laplacian, it is natural to explore the properties of these fractal Laplacians 
that are natural analogs of results that are known for the usual Laplacian. As for the 
fractal analog of the mean value property, we won't state the nature of the sets on which 
we do the averaging here, but will say that if K is a fractal set and x G K, we investigate 
whether there is a sequence of sets Bk(x) containing x with f] k B^ix) = {x} such that 



for every harmonic function u. Moreover, for general u not assumed harmonic, is there a 
formula analogous to ( 11.11) ? 

In the present paper, we will mainly deal with the Sierpinski gasket SQ. This set is 
a key example of fractals on which a well established theory of Laplacian exists [3 — 7]. 
Since the mean value property plays a very important role in the usual theory of harmonic 
functions, it is of independent interest to understand the similar property of harmonic 
functions on the Sierpinski gasket. We will prove that for each point x G SQ \ Vq, (Vq is 
the boundary of SQ.) there is a sequence of mean value neighborhoods Bk(x) depending 
only on the location of x in SQ. {Bk(x)} forms a system of neighborhoods of the point 
x satisfying f] k B k (x) = {x}. On such sequences, we get the fractal analogs of the mean 
value properties of both the harmonic functions and the general functions which belong 
to the domain of the fractal Laplacian satisfying some natural continuity assumption. We 
also investigate the extent to which our method can be applicable to other p.c.f. self- 
similar sets, but it seems that it strongly depends on the symmetric properties of both 
the geometric structure and the harmonic structure of the fractals. 

The paper is organized as follows: In Section 2 we briefly introduce some key notions 
from analysis on the Sierpinski gasket. In Section 3 and Section 4, we prove the mean 
value property for harmonic functions and general functions on SQ respectively. Section 
5 contains a further extension of the mean value property to p.c.f. self-similar fractals 
with Dihedral-3 symmetry. An interesting open question is to what extent the results of 
Section 4 can be extended to this class of fractals. See [2] for a related result concerning 
solutions of divergence form elliptic operators. 
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2 Analysis on the Sierpinski gasket 



For the convenience of the reader, we collect some key facts from analysis on SQ that 
we need to state and prove our results. These come from Kigami's theory of analysis 
on fractals, and may be found in [31 HI [5] . An elementary exposition may be found in 
[6j [7]. Recall that SQ is the attractor of the i.f.s (iterated function system) in the plane 
consisting of three homotheties {Fq, F\, F 2 } with contraction ratio 1/2 and fixed points 
equal to the three vertices {go, Qi, Q2} of an equilateral triangle. Then SQ is the unique 
nonempty compact set satisfying 



We refer to the sets F^SQ) as cells of level one, and by iterating (12. ip we obtain the 
splitting of SQ into cells of higher level. For a word w = (wi, w 2 , ■ ■ ■ , w m ) of length m, 
the set F W (SQ) = F Wl o F W2 o ■ ■ ■ o F Wm (SQ) with Wi 6 {0, 1, 2}, is called an m-cell. The 
fractal SQ can be realized as the limit of a sequence of graphs ro,ri, • • • with vertices 
Vo Q V\ C • • • . The initial graph r is just the complete graph on Vq = {qo, qi, q%}, which 
is considered the boundary of SQ. See Fig. 2.1. Note that SQ is connected, but just 
barely: there is a dense set of points J , called junction points, defined by the condition 
that x G J if and only if U \ {x} is disconnected for all sufficiently small neighborhoods 
U of x. It is easy to see that J consists of all images of {qo, qi, q 2 } under iterates of the 
i.f.s. The vertices {qo, q±, q 2 } are not junction points. All other points in SQ will be called 
generic points. In the SQ case, J = V* \ Vq, where V* = {J m V m . However, it is not true 
for general p.c.f. self-similar sets. In all that follows, we assume that SQ is equipped with 



Fig. 2.1. The first 3 graphs, To,! 1 !,!^ in the approximation to the Sierpinski gasket. 



the self-similar probability measure ji that assigns the measure 3 m to each m-cell. 
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(2.1) 



i=0 




3 



We define the unrenormalized energy of a function u on T m by 

E m (u) = ^2 ( u ( x ) - u (y)) 2 - 

The energy renormalization factor is r = |, so the renormalized graph energy on T m is 

£ m (u) = r- m E m (u), 

and we can define the fractal energy £{u) = linim^oo £ m {u). We define domS as the space 
of continuous functions with finite energy. Then 8 extends by polarization to a bilinear 
form S(u, v) which serves as an inner product in this space. 

The standard Laplacian may then be defined using the weak formulation: u G domA 
with Au = f if / is continuous, u G domS, and 



£(u, v) = — 




for all v G domoS, where dom^S = {v £ £ : v\y = 0} . There is also a pointwise formula 
(which is proven to be equivalent in [7]) which, for points in K \ Vo computes 

3 

Au(x) = - hm 5 m A m u(x), 
where A m is a discrete Laplacian associated to the graph T m , defined by 

A m u(x) = ^2 ( u (y) - u ( x )) 

for x not on the boundary. 

It is not necessary to invoke the measure to define harmonic functions, although it is 
true that these are just the solutions of Ah = 0. The more direct definition is that 

for every nonboundary point and every m. This can be viewed as a mean value property of 
h at the junction points. The space of harmonic functions is 3-dimensional and the values 
at the 3 boundary points may be freely assigned. Moreover, there is a simple efficient 
algorithm, the "| — | rule", for computing the values of a harmonic function exactly at 
all vertex points in terms of the boundary values. The harmonic functions satisfy the 
maximum principle, i.e., the maximum and minimum are attained on the boundary and 
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only on the boundary if the function is not constant. We call a continuous function h a 
piecewise harmonic spline of level m if h o F w is harmonic for all \w\ = m. 
The Laplacian satisfies the scaling property 

A(uoF0 = ^(An) o Fi 

and by iteration 

A(uoF w ) = — (Au)oF w 

for F = F o F o ■ ■ ■ o F. 

Although there is no satisfactory analogue of gradient, there is a normal derivative 
d n u(qi) defined at boundary points by 

d n u(qi) = lim V] r~ m (u(qi) - u(y)), 

m— ¥00 ' ' 

the limit existing for all u G domA. The definition may be localized to boundary points 
of cells: for each point x G V m \ Vq, there are two cells containing x as a boundary point, 
hence two normal derivatives at x. For u G domA, the normal derivatives at x satisfy 
the matching condition that their sum is zero. The matching conditions allow us to glue 
together local solutions to Au = f. 

As is shown in [31 HI E], the Dirichlet problem for the Laplacian can be solved by 
integrating against an explicitly given Green's function. Recall that the Green's function 
G(x,y) is a uniform limit of Gm{x,v) as M goes to the infinity, with Gm defined by 

M 

G M (x,y) = J2 E 9(z,z')4 m+1 \x)^ +1 \y) 

m=0 z,z>£V m+1 \V m 

and 

g(z, z) = ^r m for z G V m+1 \ V m , 

g(z, z') = ^r m for z, z' G V m+ \ \ V m with z, z' G F W (SQ) for \w\ = m, and z ^ z', 

where ifj™(x) denotes a piecewise harmonic spline of level m satisfying ifj z m \x) = S z (x) 
for x G V m . 



3 Mean value property of harmonic functions on SQ 

Lemma 3.1. (a) Let C be any cell with boundary points Po,Pi,P2> an d h any harmonic 
function. Then 

-^jM f i=^(h(p ) + h(p 1 ) + h(p 2 )). 
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(b) Let p be any junction point, and C\, Ci the two m-cells containing p. Then 

ll(r \\r\ I hdf * = 

fJ,{Ui U G 2 J J Cl uC 2 

Proof. The space of harmonic functions on C is three-dimensional. A simple basis 
{ho,hi,h2} is obtained by taking hj(pj) = 1 and hj(pk) = for k ^ j. Noticing that 
ho + hi + h2 is identically 1 on C, by symmetry, J c hidfi = \^{C) for each %. Hence 
(a) follows, (b) follows by combining (a) for C — C\ and C = Ci with the mean value 
property of h at p. □ 

Note that (b) gives a trivial solution to the problem of finding mean value neighbor- 
hoods for junction points. 




Fig. 3.1. C w and its three neighboring cells. 
The right part of the figure refers to the proof of Lemma 4.1. 

Given a point x in SQ \ V , consider any cell F w (SQ) (denote it by C w ) containing the 
point x, with boundary points F w qi = p^. Choose the cell C w small enough, such that 
it does not intersect Vq. Then it must have three neighboring cells Cq, Ci and Ci of the 
same level with Cj intersecting C w at Pi. Denote by D w the union of C w and its three 
neighbors. See Fig. 3.1. In this section, we will describe a method to find a subset B of 
D w , containing C w , such that for any harmonic function h, the mean value of h over B is 
equal to its value at x, i.e., Msih) = h(x) where Msih) is defined by 



M B (h) = 




Then we will call the set B a k level mean value neighborhood of x associated to C w where 
k is the length of w. 
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Let h be a harmonic function on SQ. The harmonic extension algorithm implies that 
there exist coefficients {aj(x)} depending only on the relative position of x and C w such 
that 

i 

Moreover, since constants are harmonic we must have 

^ ai(x) = 1 

i 

and by the maximum principle all Oj(x) > 0. Let W denote the triangle in M 3 with 
boundary points (1, 0, 0), (0, 1, 0) and (0, 0, 1) and n w the plane in M 3 containing W . So 
{(ao(x), ai(x), a,2 (x))} G W for any x G C w . However, not every point in W occurs in this 
way. 

On the other hand, given a set B such that C w cfiC D w , by linearity we have 

M B (h) = J2^h( Pi ) (3.1) 

i 

for some coefficients (00,01,02) depending only on the relative geometry of B and C w . 
Again we must have J^O; = 1 by considering h = 1. So (00,01,02) G ttw (Later we will 
show that (00,01,02) does not have to belong to W for some sets B.) Thus we have a 
map, denoted by T from the collection of B's to nw- If we can show that the image of 
the map T covers the triangle W for some reasonable class of sets B, then we can get a 
set B over which the mean value property holds for all harmonic functions. Moreover, if 
we can prove T is one-to-one, then we get a mean value neighborhood B of x associated 
to C w , that is unique within the collection of sets we are considering. 

The above is the basic idea of our method. Hence, the remaining task in this section 
is to find a suitable class B of sets B such that there is a map T from B to 7r^, such 
that T(B) covers the triangle W. Comparing with the usual mean value neighborhoods 
(they are just balls in the Euclidean case), it is reasonable to require B to be as simple as 
possible. They should be connected, possess some symmetry properties, depend only on 
the relative geometry of x and C w , and be independent of the level of C w and the location 
of C w . 

In the following, we use p to denote the distance from p to the line containing p 1 and 
P2, namely, p is the length of the height of the minimal equilateral triangle containing C w . 
Call p the size of C w . 
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Definition 3.1. Let co,ci,C2 be three real numbers satisfying < q < 1, denote by 
B(cq, Ci, c 2 ) the set 

B(co, ci, c 2 ) = C w U E U E 1 U £ 2 , 

where each Ei is a sub-triangle domain in Ci obtained by cutting Ci symmetrically with a 
line at the distance Cip away from the vertex pi. 

Remark. See Fig. 3.2 for a sketch of _B(c , c 1; c 2 ). For example, 5(0,0,0) = C w and 
5(1,1,1) = D w . Denote by 

B = {B(c , Cl ,c 2 ) :0< a < 1} 

the natural 3-parameter family of all such sets. Each member of B contains C w and is 
contained in D w . Denote by 

a:B^A 

the natural one-to-one projection with a(B(c , c\, c 2 )) = (c , ci, c 2 ) ; where A = {(c , c\, c 2 ) : 
< Ci < 1}. 




Fig. 3.2. The relative geometry of -B(co,ci,c 2 ) and C w . 

For each vector (co, ci, c 2 ) G A, there is a unique vector (ao, ai, a 2 ) G vr^ corresponding 
to the set B(co, c±, c 2 ), satisfying (13.11) where B is replaced by -B(co, ci, c 2 ). This defines a 
map T from A to n w . Then T described above from B to 7Tiy is exactly T = T o a. 

The following lemma shows that the value T(co, ci, c 2 ) is independent of the particular 
choice of C w , which benefits from the symmetric properties of the set B(cq, Ci, c 2 ). 

Lemma 3.2. T(cq,ci,c 2 ) is independent of the particular choice of C w . 
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Proof. Let h be a harmonic function. First we consider the integral j E hdfi. Denote 
by {si,ti,pi} the boundary points of Cj. By linearity, t J E hdfi can be expressed as 
a non-negative linear combination of {h(si),h(ti),h(pi)}, which by symmetry must have 
the form 

/ hdfi = (mih(pi) + ni(h(si) + hfc))) fi(C w ), (3.2) 

for some appropriate non- negative coefficients rrii,ni. Notice that in (13.21) . the coefficients 
mi,rii are independent of the location of Cj in SQ. Actually, they only depend on the 
relative position of in Cj, i.e., m i; rij depend only on q. Using the mean value property 
at pi, namely 

Ah{pi) = fifa-i) + h(p i+1 ) + h(si) + h(U), 

we obtain 

/ hdfi = {mih(pi) + ni{Ah(pi) - fifa-i) - h(p i+1 )))fi(C w ) 

JEi 

= ((Am + mi)h(pi) - ni(h(pi-\) + h(p i+ x)))fi(C w ). 

Notice that the ratio of fi(Ej) to fi(Ci) also depends only on Cj. Combined with Lemma 
3.1(a), we see that (a , ai, a 2 ) = T(c , ci, c 2 ) is independent of the particular choice of C^, 
depending only on (c , ci, c 2 ).D 

We will show the image of the map T covers the triangle W. More precisely, T(co, ci, c 2 ) 
will fill out a set PU which is a bit larger than W. Denote by P = (1, 0, 0), P\ — (0, 1, 0) 
and P 2 = (0, 0, 1) the three boundary points of the triangle W in M 3 and by O the center 
point of W. 




Fig. 3.3. a 1/6 region of W surrounded by OQo, OP2 and i-bQo- 



Lemma 3.3. T(0, 0, 1) = P 2 and T(0, 1, 1) = Q where Q = { — |, |, |} is a point in 
ii w located outside ofW. 
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Proof. From Definition 3.1, .8(0,0, 1) = C w U C 2 . Hence by Lemma 3.1(b), for any 
harmonic function h, we have Mb(o,o,i) (h) = h(p 2 ). This implies T(0, 0, 1) = P 2 . Similarly, 
5(0, 1, 1) = C w U C\ U C 2 , then for any harmonic function h, still using Lemma 3.1, we 
get 



Mb(oii)(^) — — 77n ( / hdfj,+ / hd/j — / /id/i 



15 5 

which gives T(0, 1, 1) = Q . □ 

Lemma 3.4. T({(0,c, 1) : < c < 1}) is a continuous curve lying outside of W , 
joining P 2 and Qq. (See Fig. 3.3.) 

Proof. From Lemma 3.3, by varying c continuously between and 1 we trace a 
continuous curve P 2 Qo joining P 2 and Qq. So we only need to prove the curve P 2 Qo lies 
outside of W. To prove this, we consider the set B = B(0, c, 1) for < c < 1. In this case 

B = C w U E x U C 2 . 

Given a harmonic function h, by the proof of Lemma 3.2, we have 

hdfi = ((Ani + m^hip-i) - rii(/i(p ) + h(P2)))KC w ), 



for some appropriate non- negative coefficients m\,ri\ depending only on c. 
On the other hand, we have 



hd[i = 2h(p 2 )[i(C u 

CwUC? 



by Lemma 3.1(b). 
Hence 



hdfi = / hdfi + / hd/i 

B J Ex JC W UC 2 

= (-ni/i(p ) + (4ni + mi)h(pi) + (2 - ni)h(p 2 ))n(C w ). 

The coefficient of /i(po) is always less than 0. Moreover, it equals to if and only if 
Ei = (c=0). Hence T(0, c, 1) will always lie on the outside of the triangle W as c varies 
between and 1. □ 

Now we come to the main result of this section. 
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Theorem 3.1. The map T from B to irw fills out a region W which contains the 
triangle W . 

Proof. We only need to prove the map T from B to ny/ fills out a 1/6 region surrounded 
by the line segments OQo, OP 2 and the curve P2Q0 as shown in Fig. 3.3. Then we will 
get the desired result by exploiting the symmetry. 

Consider a subfamily B\ = {-6(0, 0, c) : < c < 1} of B. If we restrict the map T 
to B\, by varying c continuously between and 1 we trace a curve (it is a line segment, 
which follows from the symmetry of E 2 ) in W joining the center O and the vertex point 
Pi- 

Consider another subfamily B 2 = {B(0, c, c) : < c < 1} of B. If we restrict the map 
T to £>2, by varying c continuously between and 1 we trace a curve (it is also a line 
segment, which follows from the symmetric effect of E\ and E 2 ) in W joining the center 
O and the point Qo across the boundary line P\P 2 with Qo located outside of W, where 
Qo is the point defined in Lemma 3.3. 

Fix a number < y < 1. Consider a subfamily C y = {B(0,c,y) : < c < y} of B. 
If we restrict the map T to C y , by varying c continuously between and y we trace a 
curve T y joining the two points T(0, 0, y) and T(0, y, y). The first endpoint T(0, 0, y) lies 
on the line segment OP 2 and the second endpoint T(0, y, y) lies on the line segment OQo- 
(See Fig. 3.3. for T y .) When y = 0, the curve r draws back to the single center point 
O. When y = 1, by Lemma 3.4, the curve Ti is a continuous curve located outside of 
the triangle W. Moreover, P 2 is the only common points of Ti and W. Hence if we vary 
y continuously between and 1, we can fill out the 1/6 region surrounded by the line 
segments OQ , OP 2 and the curve P 2 Qo- n 

Remark. In the proof of the above theorem, we actually only consider those sets B 
in B which are contained in the union of C w and subsets of only two neighbors. See Fig. 
3.4- Of course, the map T restricted to this subfamily is one-to-one, which can be easily 
seen from the proof. Hence instead ofB, the map T is one-to-one from B* onto W , where 

B* = {B(0, Cl , c 2 ) : < a < 1} U {B(c , 0, c 2 ) : < a < 1} U {B(c , c u 0) : < a < 1}. 

Based on the discussion in the beginning of this section, we then have 

Theorem 3.2. For each point x G SQ \ V , there exists a system of mean value 

neighborhoods Bk(x) with f) k Bk(x) = {x}. 

Proof. Let ko be the smallest value of k such that there exists a k level cell C w 

containing x but not intersecting V . (k depends on the location of x in SQ.) Then 
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Fig. 3.4. The 3 shapes of B E B* associated to C w shown in Fig. 3.1. 

by using Theorem 3.1 we can find a sequence of words w^ k ' of length k (k > k ) and a 
sequence of mean value neighborhoods B k (x) associated to C w (k). Obviously, {B k (x)} k > ko 
will form a system of neighborhoods of the point x satisfying f] k>ko B k (x) = {x}. □ 

4 Mean value property of general functions on SQ 

In this section, we extend the mean value property to more general functions on SQ. 
Given a point x in SQ \ Vq and a cell C w = F W (SQ) containing x, for each mean value 
neighborhood B of x associated to C w , we assign a constant c B to B. We want 

M B {u) — u(x) « c B Au(x) 

for u in domA. More precisely, let {B k (x)} k > ko be the system of mean value neighborhoods 
of the point x; we want 

l im -J— {M Bk(x) - u{x)) = Au(x) (4.1) 

fc^oo C Bk ( x ) 

for appropriate functions in the domain of A, which is the desired fractal analog of (II. ip . 

For this purpose, let v be a function on SQ satisfying Av = 1. For each point x in 
SQ \ Vq, and each mean value neighborhood B of x, define c B by 

c B = M B (v) - v(x). 

Note that the result is independent of which v, because any two such functions differ by 
a harmonic function and the equality M B (h) — h(x) = always holds for any harmonic 
function h. So we can choose 




G(x,y)dfi(y), 
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which vanishes on the boundary of SQ. Here G is Green's function. 

We will prove that eg is controlled by the size of C w . More precisely, we will prove: 
Theorem 4.1. Let x G SQ \ V and B be a k level mean value neighborhood of x. 

Then 

1 1 

c ^ <c B < Cl -£ 

for some constant Cq,Ci which are independent of x. 

To prove Theorem 4.1, we need the explicit expression for the function v. Recall from 
Section 2 that v(x) is the uniform limit of %(x) for 



v M (x) = - J G M (x,y)d/j,(y). 

Interchanging the integral and summation, 

m „ 
Vm{x) = - E E 9(z, z') I ^ +1 \yW(y)^ m+1 \x). 

m=0 2 , 2 'eV m+1 \V m 

Notice that for each z G V m+ i \ V m , ifj z m+1 ^ is a piecewise harmonic spline of level 
(m + 1) satisfying ip z m+1 \y) = 5 z (y) for y G V m +i- More precisely, ipi m+1 ^ is supported in 
the two (m + l)-cells meeting at z. If F T (SQ) is one of these cells with vertices z, z\ and 
z 2 , then ijji m+1 ^ + ipi™ + restricted to F T (SQ) is identically 1. Thus 



JF T (sg) 6 



By symmetry all three summands have the same integral, so j F , S g, dfi = 3^+2. 

Together with the contribution from the other (m + l)-cell we find for each z G V m+ \ \ V m , 

^ +1) (yW(y) = 3^5. (4-2) 

Hence 

2 M 1 

m=0 z,z'£V m+1 \V m 

Substituting the exact value of g(z, z')(see Section 2 and details in [7] page 50) into it, we 
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get 



vmIx 



\Eh E E ^.^< m+ "M 

m=0 \| W |=m 2 , 2 'eF„(y )\F„(yi) ; 

&(e e (Jr+2|r)^>(,)) 

1 M 1 



15 ^5 

m=0 



for 



Thus 



<f>m(x) = ^ ^i"" 



zGV m+ i\y„ 



15 ^ 5 r 

m=0 

Remark. T/ie function v is invariant under Dihedral-3 symmetry. 

This is a direct corollary of the fact that each 4> m (x) is invariant under D 3 symmetry. 

Due to the above remark, we may assume that D w associated to C w has a fixed shape 
as shown in Fig. 3.1 without loss of generality. We now show that although c B depends 
on the relative position of x in C w , it does not depend on the location of x or C w in SQ. 

Lemma 4.1. Let x,x' be two distinct points in SQ \ Vq. Let C w and C w < be two k 
and k' level neighboring cells of x and x' respectively. Denote by B and B' two mean 
value neighborhoods of x and x' respectively. If B and B' have the same shapes (the same 
relative locations associated to C w and C w > respectively), then 

c B = 5 k '- k c B ,. 

In particular, if B and B' have the same levels and same shapes, then cb = cb> ■ 

Proof. D w can be decomposed into a union of a k level cell D$ and a [k — 1) level 

cell L>i 2) as shown in Fig. 3.1. Denote by q the junction point connecting D$ and D$ ■ 

Similarly, D w > can also be written as a union of a k! cell £>$ and a (k' - 1) cell D$ with 

a junction point q' connecting them. 

Let r be the linear function mapping D w onto D w /. Suppose Dw = F a (SQ) and 

D$ = Fp(SQ) where a and f3 are the corresponding words of and D$ respectively. 

Similarly, denote by a' and (3' the corresponding words of D^) and D^). Hence we can 
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write r as r(z) = F a > o F'^z) if z G and r(z) = Fp o if 2 6 £>i 2) . In 

particular, r(g) = </ and r(x) = x' . 

Consider the function (v o F a — 5 k '~ k v o F a >) defined on SQ. Noting that \a\ = k and 
\a'\ = k', using the scaling property of A(see details in [7] page 33), we have 

A(v o F a — 5 k '~ k v o F a >) = r |a| ^At) o F a - h k '^ k r\ a '\^Av o F a > = 0, 

which shows that the difference between v o F a and 5 k '~ k v o F a > is a harmonic function. 
Hence the difference between v and 5 k '~ k v o r on D$ is harmonic. A similar discussion 
will show that the difference between v and 5 k '~ k v or on is also harmonic. Since the 
matching condition on normal derivatives of (v — 5 k '~ k v o r) at q holds obviously, we have 
proved that A(v — 5 k '~ k v o r) = on D w , i.e., the function {v — 5 k '~ k v o r) is harmonic 
on D w . 

By the definition cb = Mb{v) — v{x) and cb> = Mb'(v) — v(x'). Notice that for the 
second equality, by changing variables we can write cb> = Mb{v o t) — v o t(x). Hence 

c B - 5 k '- k c B i = M B {v - h k '~ k v or) - (v - 5 k '~ k v o r)(x) = 0, 

since (v — 5 k '~ k v o r) is a harmonic function on D w . □ 



Proof of Theorem 4-1- 
Estimate of cb from above. 

From Lemma 4.1, since cb depends only on the relative geometry of B and C w and 
the size of C w , but not on the location of C w , we may assume that D w is contained in a 
{k — 2) level cell C in SQ without loss of generality. 

By the definition of Cb, we may write 

c B = M B (v) - v(x) = lim ( — — - / v M dfj, - v M (x)) . 

M^oo \fl(B) J B J 

Substituting the exact formula of % into it, we get 

I 00 I 

m=0 

for 



ZdVm+l\Vm 

Notice that each <p m is a piecewise harmonic spline of level m+1. So when m+1 < k—2, 
is harmonic in the cell C, which yields that M B (<f) m ) — <f> m (x) = 0. So the first k — 2 



15 



terms in the infinite series of v will contribute to cb- Hence 

1 00 1 

It is easy to see that this implies 

1 11/" 

|c fl | < — ^ ^^b) y B 1^(2/) " Mx)W{y). 



m=k—2 

Then by the maximum principle, we finally get 



1^1 25 1 

m=k— 2 



Estimate of eg from below. 




Pi P2 

Fig. 4.1. a 1/3 region of C w . 

Without loss of generality, we assume that x is located in the 1/3 region of C w as 
shown in Fig. 4.1, i.e., x is contained in the triangle T PltP2)0 , where o is the geometric 
center of C w . Then by the proof of Theorem 3.1, B is a subset of the union of C w and two 
of its neighbors C\ and CV Hence we can write B = C w U E\ U E2, where E{ = B n C{. 

Claim 1. Let B = F (SQ)UEiUE2, where Ei is a triangle obtained by cutting F^SQ) 
symmetrically with a line below the top vertex Fiq .(see Fig. 4- If B and B have the 
same shapes, then 



c B = 5 1_fc c 



B- 



This is a direct corollary of Lemma 4.1. 

We only need to prove that eg for B defined in Claim 1 has a positive lower bound. 
For simplicity of notation, in all that follows, we write B instead of B. In other words, 
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we only need to consider B whose associate cell C w is F (SQ). In this setting, pi = F qi, 

C t = Fi(SG) and c 2 = F 2 (sg). 

We write v = — v where v is the non- negative function defined by 

m=0 

For each M > 0, denote by 

the partial sum of the first M + 1 terms of t>. Then % converges to t> uniformly as 
M ->■ oo. 

We have the following three claims on 7>. 

Claim 2. < v < 1 on and n ia/ces constant 1 a/ong t/ie maximal inner upside- 
down triangle contained in SQ . 

Proof. Consider the partial sum function %. Obviously, % is a (M+l)-level piecewise 
harmonic function on 5^. For convenience, denote by V the maximal inner upside-down 
triangle contained in SQ. We divide the vertices Vm+i into three parts, V^ +1 , V^ +1 and 
Vm+u where V' M+ i consists of those vertices lying along V, V^ +l consists of those vertices 
at distance 2~( A/+1 ) from V, and V^ +l consists of the remain vertices. Then by using the 
"| — |" rule, an inductive argument shows that Vm = 1 on V' M+1 , % = 1 — -^kt on V^ +1 , 
and Vm < 1 — 517 on V^ +1 . Since v is the uniform limit of % and V^ //+1 goes to V as M 
goes to the infinity, we then have < v < 1 on SQ and v = 1 on V. □ 

Claim 3. For each x contained in the triangle T pi P2 , v(x) > ||. 

Proof. For r = (0, 1, 1), (0, 1, 2), (0, 2, 1) and (0, 2, 2), by using the harmonic extension 
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algorithm, namely, the "| — |" rule, we get that 



2 

v(F T q ) = v 2 (F T q ) = £ ^0 m (F r g o ) = 1.1 + 1| + 1.1 = ?1 



m=0 

where |, | and 1 are the values of O , 0i and 4> 2 at F r g respectively. Also, for those r, 
by Claim 2, we have 

u(F T gi) = v 2 {F T qi) = v(F T q 2 ) = v 2 (F T q 2 ) = 1. 

Notice that for each point x in the triangle T Pl)P2)0 , x is contained in one of the four 3- 
level cells F n(SQ), F 012 (SQ), F 021 (SQ) and F 022 (SQ). Since v 2 is harmonic in each such 
cell, by using the maximal principle, we get that 

24 

v 2 (x) > — . 
v 1 ~ 25 



Hence v(x) > || since each term in the infinite series of v is non-negative. □ 
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Claim 4. Mb(v) < g. 
Proof. First of all we prove that 



/ v(y)dfM(y) = 1. 



'Fo(SC) 

We need to compute J Fo t S g^ 4>m{y)d^{y) for each non-negative integer m. For each m > 0, 

<f>m(yW(y) = \ ■ s m+1 • = I 

F (sg) 6 6 y 

by using (14. 2 p and the fact that m = ^2 zeVm+1 \ Vm ^z +V> • Hence 



/" 2 V °° A 1 5 

/ WM?/) = qE^ = Tl 



v{y)dfi{y) = ~2^- 

By our assumption, the mean value neighborhood B can be written as 

B = F (Sg)U E 1 \J E 2 , 
where E* — B D C, . Hence we have 



Mb(v) = I v(y)dfi{y)+ I v{y)dfi(y) + / v(y)d/i(y) 

>F (Sg) J E x JE 2 



KB) 



< "Vm [ I v(y)dfi(y) + I 1 • dn{y) + / 1 • t^(j/) 



5/18 + + fi{E 2 ) 

1/3 + + fi(E 2 ) ' 
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where the inequality follows from Claim 2. Since < /i(£'i)+/i(£ , 2 ) < ~, ■Jf^ is increasing 

3 + X 

in x ^ 

5/18 + fijEi) + /i(E 2 ) 5/18 + 2/3 _ 17 
l/3 + /i(£i)+//(£; 2 ) - 1/3 + 2/3 ~ 18' 

Hence we always have 

M B (tO<— . □ 
Now we turn to estimate c B . Obviously, 

c B = M B (v) - v(x) = ~ (Mb(v) - v(x)) . 

15 

By Claim 3 and Claim 4, we notice that M B (v) — v(x) < ^ — || = — Hence 

1 7 

c > > o. □ 

15 450 

On the other hand, given a point x and C w = F W (SQ) a fc level neighborhood of x, for 
any u G domA, we write 

it = + (Au{x))v + i? (fc) 
on C^, where hP^ is a harmonic function defined by 

h {k) + (Au(x))v\ dCw = u\ dCw . 
It is not hard to prove the following estimate: 

Lemma 4.2. Let u G domA with g = Au satisfying the following Holder condition 

\g{y)-g{x)\<n k , (0 < 7 < 1) 
for all y G C w . Then the remainder satisfies 



i? (fc) =o(( 



on C w (hence also on Bk(x)). 

Proof. It is easy to check that AR^(y) = Au(y) — Au(x) and PS k \y) vanishes on the 
boundary of C w . Hence RS k ^ is given by the integral of Au(y) — Au(x) on C w against a 
scaled Green's function. Noticing that the scaling factor is (\) k and 



\Au(y) - Au(x)\ < ey k 



we then get the desired result. □ 
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This looks like a Taylor expansion remainder estimate of u at x. See more details on 
this topic in [8]. 

Remark. // we require u G domA 2 , then the remainder satisfies 

,3 1, 



i? (fc) = O 



~) k 

5 5 ; 



on C w (hence also on Bf.(x)). The reason is that in this case Au satisfies the Holder 
condition that \Au(y)—Au(x)\ < c(|) fc for ally G C w , because A 2 u is assumed continuous, 
see 0, Theorem 8.4- 

Using the above lemma and Theorem 4.1, we then have the following main result of 
this section. 

Theorem 4.2. Let u G domA with g = Au satisfying the Holder condition \g(y) — 
g(x)\ < C7 fc for some 7 with < 7 < 1, for all x,y belonging to the same k level cell. 
Then 

lim (M Bk (x){u) - u(x)) = Au(x). 

k^oo C Bk ( x ) 

Proof. Using Taylor expansion of u and noticing that MB k ( x )(h^ k ' > ) — hS k \x) = 0, 
M Bk ( x )(v) - v(x) = c Bk ( x ), we have 

— (M Bk(x) (u) - u(x)) - Au(x) = — (M Bk(x) (RW) - RW(x)) 

C B k (x) C Bk(x ) 

= —o(a) k ) = o( 1 k ). 

CB k (x) V 5 / 

Hence letting k — > 00, we get the desired result. □ 



5 p.c.f. fractals with Dihedral-3 symmetry 

The results for SQ should extend to other p.c.f. fractals which possess symmetric prop- 
erties of both the geometric structure and the harmonic structure. We assume that a 
regular harmonic structure is given on a p.c.f. self-similar fractal K. The reader is re- 
ferred to [U [7j for exact definitions and any unexplained notations. We assume now that 
JjVo = 3 and all structures possess full D 3 symmetry. This means there exists a group 
Q of homeomorphisms of K isomorphic to D% that acts as permutations on Vo, and Q 
preserves the harmonic structures and the self-similar measure. 

Assume that the fractal K is the invariant set of a finite iterated function system of 
contractive similarities. We denote these maps {-Fi}i=i,... ,tv with N > 3. Let denote 
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the 2-th resistance renormalization factor and p,i denote the i-ih weight of the self-similar 
measure /x on K. In general, it is not necessary that all r^'s and all /ij's be the same, but 
here we must have ro = T\ — • ■ ■ — rjy and /Jo = A*i = • • • = /-tjv from the above Dihedral-3 
symmetry assumption. We denote Vq = {Qo, Qi, Q2} the set of boundary points. 

Examples, (i) The Sierpinski gasket SQ. In this case all = 3/5 and all \ii = 1/3. 

(ii) The hexagasket, or fractal Star of David, can be generated by 6 maps with simul- 
taneously rotate and contract by a factor of 1/3 in the plane. Thus Vq consists of 3 points 
of an equilateral triangle, and V\ consists of the vertices of the Star of David, as shown in 
Fig. 5.1. Although the same geometric fractal can be constructed by using contractions 
which do not rotate, this gives rise to a different self-similar structure (in particular with 
JjVo = 6). Our choice of self-similar structure destroys the D$ symmetry of the geometric 
fractal, but it has the advantage of easier computation. In this case, all rj = 3/7 and all 
Hi = 1/6. Note that in this example there exist points in V\ that are not junction points, 
(iii) The level 3 Sierpinski gasket SQ3, obtained by taking 6 contractions of ratio 1/3 as 




r ii 

Fig. 5.1. The first 2 graphs, ro,Ti in the approximation to the hexagasket. 

shown in Fig. 5.2. Here we have all r, = 7/15 and /i, = 1/6. Note that all seven vertices 
in V\ \ Vo are junction points, but the one in the middle intersects three 1-cells. In a 
similar manner we could define SQ n for any value of n > 2. 

We prove that there are results analogous to Theorem 3.1, which yield the existence 
of mean value neighborhoods associated to K. 

Given a point x in K\Vo, consider any cell F W K = C w with boundary points po,Pi,P2 
containing the point x. Without losing of generality, we may require that the cell C w 
does not intersect V . For each i, denote by C^i, ■ • • , the neighboring cells of C w of 
the same size, intersecting C w at pi, where U is the number of such cells. It is possible 
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Fig. 5.2. The graph of the V\ vertices of the level 3 Sierpinski gasket. 



that li = for some i since Pi may be a non-junction point. If this is true, the matching 
condition says that the normal derivative of any harmonic function h must be zero at this 
point, which yields that the value of h at this point is the mean value of the values of 
h at the other two boundary points of C w . In other words, the restriction of all global 
harmonic functions in C w is two dimensional. Denote by D w the union of C w and all its 
neighboring cells, i.e., 



Two cells C w and C w i are said to have the same neighborhood type if they have the same 
relative geometry with respect to D w and D w i respectively. It is obvious that there only 
exist finitely many distinct types. For example, for SQ, all cells have exactly only one 
neighborhood type. For SQs, the number of the finite types is 3. For SQ n (n > 4), the 
number of the finite types becomes 4. For the hexagasket gasket, the number of the finite 
types is 2. 

Let h be a harmonic function on K. Given a set B containing C w , define 



for some coefficients (00,01,02) satisfying = 1. Notice that this is true for SQ. In 
that setting, a harmonic function is uniquely determined by its values on the boundary of 



D w — C w U Cjj. 




the mean value of h over B. We are interested in an identity 




(5.1) 
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any given cell C w because the harmonic extension matrix associated with C w is invertible. 
However, in the general case, the harmonic extension matrices may not be invertible. So 
we can not prove (15.11) for every set B simply by linearity. However, it will suffice to show 
that the equality (15.1 ft holds for certain specified sets B. 

Consider a set B which is a subset of D w , containing C w . Then B must be made up 
of four parts, i.e., 

B = C w U E U E 1 U E 2 

where Ei = B fl C, with = U l j =1 C it j. It is possible that Cj may be empty since Pi 
may be a nonjunction point. We can also subdivide each Ei into U small pieces, i.e., 
Ei = UjEij for Eij = Ei fl Cij. For each i, we require that E^i, • • • , be of the 
same size and shape. Moreover, in analogy with the SQ case, we require that each Eij to 
be a symmetric (under the reflection symmetry that fixes pi) cutoff sub-triangle of Cij, 
containing pi as one of its vertex points. This means that there is a straight line Ljj, 
symmetric under the reflection symmetry fixing pi, cutting C%j into two parts, and E^j 
is the one containing pi. For each Eij, define the distance between pi and the line Ljj 
the size of Eij. Of course, for each fixed i, E it i, ■ ■ • , E^ have the same sizes. We call the 
common value the size of E^ Suppose the size of every Cjj is p. (Of course, they are all 
equal.) Then for each i, the size of Ei is Qp where the coefficient < q < 1. Hence we 
can write the set B = B(cq, c±, C2). (If p% is a nonjunction point, then q should always be 
0.) For example, suppose that the boundary points of C w consist of junction points, then 
B(0, 0, 0) = C w and 5(1, 1, 1) = D w . Denote by 

B = {B(c , Cl , c 2 ) : < a < 1} 

the family of all such sets. Then we can show that the formula (15. ip holds for each B £ B. 

Proposition 5.1. Let B e B, then for any harmonic function h, we have (15. ip for 
some coefficients (00,01,02) independent of h. Moreover, Yli a i = 1- 

Proof. Each B G B can be written as B = C w U E U E\ U Given a harmonic 
function h on K, for fixed z, we first consider the integral J E hdfi. Obviously, 

/ hd/j, = / hd\i. 

J Ei j J Eij 

For each 1 < j < l iy denote by {z i: j, Wij,pi} the boundary points of C it j. Since each is 
contained in Cij, hdp: can be expressed as a linear combination of o(pj), h(zij) 

and h(wij) with non- negative coefficients independent of the harmonic function h. Since 
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the set Eij is symmetric under the reflection symmetry fixing pi, the two coefficients with 
respect to h(zij) and h(wij) must be equal. In other words, we can write 

/ hd/j, = (m itj h(pi) + n itj h{z itj ) + n itj h(w itj )) fi(C w ) 

for rriij,nij > 0. Moreover, since for each fixed i, Eij are in the same relative position 
associated to Cij for different j's, j E hdfi can be expressed as a linear combination of 
h(pi), h(zij), h(wij) with the same coefficients for different j's. Hence we can write 

J hdfi = ^rriih^pi) + n» ^ (h( z i,j) + K w i,j))j K C w), 

for suitable coefficients mi, Hi > 0. The mean value property at the point pi says that 

+ h{w id )) = (2k + 2)h(pi) - (/ife-i) + 

i 

Combining the above two equalities, we get 

hdfi = ((mi + 2kni + 2rii)h(pi) - n^p^i) - nih(p i+1 )) fi(C w ). 

On the other hand, by the linearities and symmetries of both the harmonic structure 
and the self-similar measure, 

hdn = (h( Po ) + h( Pl ) + h(p 2 )) . 



Since the ratio of (i(Eij) to fi(C w ) depends only on q, we have proved that Ms(h) 
can be viewed as a linear combination of the values of h on the boundary points of C w , 
i.e., 

M B (h) = Y, a * h &), 

i 

where the combination coefficients are independent of h. Moreover, we must have = 1 
by considering h = 1. □ 

Remark 1. This means that Ms(h) is a weighted average of the values h(po),h(p\) 
and h(p 2 ). Moreover, if one of the boundary points, for example p 2 , is a nonjunction 
point, then by the fact that h(p 2 ) = \(h(po) + h(pi)), we have 

M B (h) = a h(p ) + axhipx) + ^a 2 (h(p ) + h(p x )) = a h(p ) + a 1 h(p 1 ) 
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for clq = do + \a,2 and oi — a\ + |«2- We o/so oai>e ao + ai = 1. Hence in this case, we 
can also view Ms{h) as a weighted average of the values of h(p ) and h(pi). 

Remark 2. The proof of Proposition 5.1 shows that (00,01,02) depends only on the 
neighborhood type of C w and the relative position of B associated to C W! and does not 
depend on the particular choice of C w . In other words, if we consider a cell C w with a 
given neighborhood type, then for each set B £ B with the expression B = B{cq, 0^,02), 
the coefficients (00,01,02) of B depend only on (cq, 01,02). 

The following is the main result in this section. 

Theorem 5.1. Given a point x £ K\V$, let C w be a cell containing x, not intersecting 
Vq, and let D w be the union of C w and its neighboring cells of the same size. Then 
there exists a mean value neighborhood B of x satisfying C w C B C D w . Moreover, for 
each point x £ K\Vq, there exists a system of mean value neighborhoods B k {x) with 
f] k B k (x) = {x}. 

Proof. We need to classify the distinct neighborhood types into three cases according 
to the number of no nj unction points in the set of boundary points of C w . 
Case 1. All boundary points of C w are junction points. 

This case is similar to what we have described in the SQ setting. Let W denote the 
triangle in R 3 with boundary points P = (1,0,0), Pi = (0,1,0) and P 2 = (0,0,1) and 
iiw the plane containing W. Notice that from Proposition 5.1, (oo,ai,a 2 ) £ irw for each 
B. We use T to denote the map from B to ttw- From Remark 2 of Proposition 5.1, the 
map T is uniquely determined by the neighborhood type of C w . Let B* be a subfamily 
contained in B defined by 

B* = {5(0, ci, ca) : < a < 1} U {P(c , 0, c 2 ) : < a < 1} U {P(c , a, 0) : < a < 1}, 

i.e., those elements B in B which have the decomposition form B = C w U Pi U E 2 or 
B = C w U P U P 2 , or B = C w U P U E 1 . Then we have 

Claim 1. The map T from B to nyy fills out a region W which contains the triangle 
W . Moreover, T is one-to-one from B* onto W. 

Proof. The proof is similar to the SQ case. The only difference is the line segments 
OQo and OP2 described in the proof of Theorem 3.1 may become continuous curves OQo 
and OP2 in the general setting. □ 

Case 2. There is one nonjunction point (for example, p 2 ) among the bound- 
ary points of C w . 

In this case, there is no neighboring cell intersecting C w at the point p 2 - Hence P 2 will 
always be empty. So B = {P(cq, C\, 0) : < q < 1} for this case. 
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As shown in Remark 1 of Proposition 5.1, for any harmonic function h on K, B G B, 
M B {h) is a weighted average of h(p ) and h(pi), i.e., 

M B (h) = a h(p ) + a-Jiipx) 

with a , a\ independent of h, satisfying a + «i = 1. Let / denote the line segment in 
M 2 with endpoints Pq = (1, 0), Pi = (0, 1) and pi the line containing /. Notice that from 
Remark 1 of Proposition 5.1, (ao, a±) G pi for each B. We still use T to denote the map 
from B to pi. From Remark 2 of Proposition 5.1, the map T is uniquely determined 
by the neighborhood type of C w . We may write T(B(c , ci, 0)) = (a ,ai) for each set 
-B(co, Ci, 0). We will show the image of the map T covers the line segment /. Similar to 
Case 1, let B* be a subfamily contained in B defined by 

B* = {B(c , 0, 0) : < c < 1} U {B(0, c x , 0) : < c x < 1}, 

i.e., those elements B in B which have the decomposition form B = C w U E or B = 
C w U Ei . Then we have 

Claim 2. The map T from B to pi fills out the line segment I . Moreover, T is a 
one-to-one map on B* . 

Proof. The proof is similar to Case 1. Denote by O = (|, o) the midpoint of /. We 
only prove the map T from B to pi fills out half of the line segment /. Then we will get 
the desired result by symmetry. 

Let hbe a harmonic function on K. We consider T({(B(c, 0, 0)) : < c < 1}). When 
c = 0, B(0, 0, 0) = C w and M Cw (h) = |(/i(po) + %i) + Kp 2 )). Combining this with the 
fact that 

h(p2) = -(% ) + %i)), 

we get 

M Cw (h) = ±(h( Po ) + h( Pl )). 

Hence T(B(0,0,0)) is the midpoint O of /. When c = 1, 5(1,0,0) = C w U C , and an 
easy calculation gives that M Cw uc — h(po). Hence T(5(l,0,0)) is the endpoint P . So 
if we vary c continuously between and 1, we can fill out the line segment joining O and 
Pq, which is half of /. □ 

Case 3. There are two nonjunction points (for example, pi and p 2 ) among 
the boundary points of C w . 

In this case, let h be any harmonic function on K. By the matching condition on both 
points pi and p 2 , h must be constant on the whole cell C w . Hence for every point x G C w , 
we could view C w itself as the mean value neighborhood of x. 
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Hence the proof of Theorem 5.1 is completed by using a same argument as that of 
Theorem 3.2. □ 

We should mention here that the result can also be extended to some other p.c.f. frac- 
tals including the 3-dimensional Sierpinski gasket. However, it seems that some strong 
symmetric conditions of both the geometric and the harmonic structures should be re- 
quired. 
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